February 14, 2001 



The Baxter Equation for Quantum Discrete Boussinesq Equation 



Kazuhiro Hikami f\ 

Department of Physics, Graduate School of Science, 
University of Tokyo, 
Hongo 7-3-1, Bunkyo, Tokyo 113-0033, Japan. 

(Received: ) 
ABSTRACT 

Studied is the Baxter equation for the quantum discrete Boussinesq equation. We explicitly 
construct the Baxter Q operator from a generating function of the local integrals of motion 
of the affine Toda lattice field theory, and show that it solves the third order operator- valued 
difference equation. 
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1 Introduction 

In analysis of the integrable models, realized is the importance of the Baxter equation since it was 
applied to the eight-vertex model (XYZ model) to compute the spectrum [[I]]. The Baxter equation 
is the functional difference relation among the transfer matrix t(X) and the auxiliary function 2(A) 
which is called the Baxter Q operator, and is powerful enough to determine the spectrum of the 
transfer matrix only from certain analytic conditions of these functions. 

Later Gaudin and Pasquier explicitly constructed the Baxter Q operator in an integral form for the 
quantum periodic Toda lattice [Q]. It was also shown [Q,|3|] that in the classical limit the logarithm 
of the Q operator is the generating function of the Backlund transformation, i.e., the Backlund trans- 
formation is given by the classical limit of the evolution equation, O — ► Q O Q^ 1 (see, for review, 
Ref. ^|). This correspondence suggests that the integrable dynamical system on space-time lattice can 
be constructed by use of the Q operator []|, §J. Recent studies indicate that the Q operator is useful 
enough to construct the eigenf unction of the Toda lattice [Q]. 
^E-mail: hikamiSphys . s . u-tokyo .ac.jp 
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Since then the explicit form of the Baxter Q operator has been derived for other integrable models 
such as the XXX model [1810], the dimer self-trapping model [10], the quantum KdV model [11], and 



the Volterra model (discrete KdV equation) [12, 13]. All these models are governed by the quantum 
algebra associated to s^2- In these cases as was pointed out in studies of the integrability of the 
conformal field theory [|Tl|], the Baxter equation may be related with the universal TZ matrix, TZ 6 
U q (si2) U q (sl2), with some quantum finite-dimensional representation for the first U q {sl2) algebra 
and an auxiliary m/rote-dimensional g-oscillator representation for the second part. In fact proposed 



was the universal procedure [14] to derive the Baxter equation from the universal TZ matrix by use of 



the g-oscillator representation. Following this strategy it was shown that the Q operator [)12| for the 
Volterra model can be constructed from the universal TZ matrix for the quantum algebra U q (si2) with 



the g-oscillator representation (15]. 

In this paper we construct the Baxter Q operator for the quantum discrete Boussinesq equation 
which has U q (slz) symmetry. In the classical limit the discrete Boussinesq equation is given as 

= -W n + W n -i + L n [L n+1 - L n _i), 

(1.1) 

= W n (L n+ 2 — L n -x). 



dt 
dW n 



dt 



We stress that L n and W n do not mean the Fourier modes of the generators L(x) and W(x) but the 
discretization thereof. We note that the higher flow of eq. ( |j~l| ) reduces in a continuum limit [ |i"6| ] to 



du 
dt 
dw 
~dt 



2 2 

-w xx + -u xxx + -u Ul 



(1.2) 



This set of equations gives the Boussinesq equation which is generated from the pseudo-differential 
operator, L = <9 3 + u d + w. In this sense we call eq. ( |1.1| ) as the discrete Boussinesq equation. 
The discrete Boussinesq equation ( |l.l[ ) was first introduced in Ref. 17 from studies of a discretization 
of the W3 algebra, motivated from the close connection between the soliton equations and the con- 
formal field theory JT8| — ^0|] - The integrability of the discrete Boussinesq equation follows from the 
fundamental commutation relation, 

{T(x) f T(y)} = [r(x, y) , T(x) (g) T(y)], (1.3) 

where we omit the explicit form of the Z3 symmetric classical r matrix. The monodromy matrix T(x) 



is given from the Lax matrix as Q16| ] 

T(x)=n^(*), 



/ x 2 -x L n+1 W n \ 



1 





1 





J 



(1.4) 
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As a result we can construct the integrals of motion from the transfer matrix Tr T(x); some of them 
are 

n 

Hi = L n , 

n 

H 2 = Y, (I i L n) 2 + L n L n+l - W n ) . 

n 

As is well known, the discrete Boussinesq equation has also an intimate relationship with a lattice 
analogue of the affine Toda field theory. We have constructed a generating function Q of the local 
integrals of motion of the affine Toda lattice field theory and have also shown that this also 
commute with the transfer matrix of the discrete Boussinesq equation. The function Q is constructed 
in terms of the quantum dilogarithm function as the intertwining operator of the vertex operators. 
As was conjectured there, we shall show here that this generating function is indeed the Baxter Q 
operator for the discrete Boussinesq equation. 

This paper is organized as follows. In Section || we briefly review the integrable structure of the 
quantum discrete Boussinesq equation. Recalling the relationship with the affine Toda field theory, 
we construct the transfer matrix and the generating function Q(A) of the local integrals of motion 
for the affine Toda lattice field theory. In Section ^ we explicitly show that the generating function 
Q plays a role of the Baxter Q operator for the transfer matrix of the discrete Boussinesq equation. 
Based on the duality of the affine Toda lattice field theory, we show in Section |] that we have the 
dual Baxter equation for the discrete Boussinesq equation. In Section || we briefly comment on the 
continuum limit of the discrete Boussinesq equation, and see a relationship between the previously 
known results. In Section ^ we consider the generalization for the A -reduced discrete KP equation 
which has the quantum group U q (s£]\[) structure. Last section is devoted for the concluding remarks 
and discussions. 

2 Quantum Boussinesq Equation 

2.1 Lattice Vertex Operator and Transfer Matrix 

To define the quantization of the discrete Boussinesq hierarchy (|L~4|), we use two sets of the discretized 
free fields 4>n ^ and (jffl , which respectively correspond to the Toda fields associated to the simple roots 
a.\ and a.2- We use the free field associated to the maximal root as 

^ 0) = -^-4 2) - (2-D 
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The commutation relations are denned as a discrete analogue of the usual relations of the free chiral 
fields by[] 

[(/>$,</>$]= 2i % form >n, a = 1,2, 



17, for n < m, 
—17, for n > m. 

Throughout this paper we set the deformation parameter q as 



(2.2) 



(2-3) 



(2.4) 



Hereafter for our convention we often use a difference of the lattice free field 

whose commutation relations become local. 

With these settings, we define the discrete screening charges as in the case of the conformal field 



theory (see e.g. Ref. 23); 



n 

for a = 0, 1, 2. We can see that the Serre relation of the algebra U q (s£s) is fulfilled; 

( Q W) 2 Q(b) _ ( ? + . Q(a) Q(6) Q(a) + q(6) (q(«))2 = Qj 



(2.5) 



(2.6) 



for a ^ b. The generators of the quantum W3 algebra are then defined so as to commute with these 
screening charges. Note that the dual screening charges as a summation of the dual vertex operators, 



:(a) 



n 

commute with Q^ a \ and that they satisfy the same Serre relation (|Z6|) replacing g with 



g = e 



(2.7) 



(2.8) 



As was shown in Ref. 16, the Lax matrix (1.4) for the discrete Boussinesq equation is gauge- 
equivalent with 

/ e iA0W-|A^ \ / e A e~ A 1 \ 

MA) = (7(A) e iAtf) + f A0£ 2 > 

e 3' 



V 



S(A) • C n • X(A), 



^A0W-lA^ 



1 e A e~ A 



1 



/ 



(2.9) 



We note that the dynamical variables a n and f3 n in Ref. ^ are defined by 



/3 n = e 



with q = e 17 . 



where we have modified the spectral parameter by = 1 — e~ 3A from eq. (1.4). The function 
g(X) is a normalization function which will be defined below. We introduce diagonal elements of the 
matrix C n as 

C n = diag(c«,cl 2) ,<i 3) 



(2.10) 

which corresponds to the weights in the vector representation 3 of si^. Following the standard way 
of studies in the integrable systems, we define the monodromy matrix and the transfer matrix as 



T(A) = ;QMA), 



ti(A) = TYT(A). 



(2.11) 



The monodromy matrix satisfies the Yang-Baxter equation [22], 



R(A - y) ■ T(A) T(/i) = T(/i) T(A) • R(A - M ). 



(2.12) 



Here the R-matrix is Z3 <g> Z3 symmetric, and reads as 



with 





( 

b 

b 


c 


c 


R(A) = 


c 


b 

a 

b 


c 




c 

{ 


c 


b 

b 

a 



\ 



a = sh(- A 17), 

6 = e^sh(|A), 



_i\ . , 1 . . 
c = e 2 sn( — 17), 



6 = e-^sh(^A), 



1 \ , , 1 . . 
c = ea sh(— — 17). 



It follows that the transfer matrix commute with each other, 

MX) , t 1 ( J j.)] = o, 

which supports the quantum integrability of our system. 

We note that under the quasi-periodic boundary condition, 



(2.13) 



(2.14) 



J n+L 



b^+L-iP^, for a =1,2, 



(2.15) 
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the transfer matrix (2.11) can be rewritten as 



L 



ti(A) = Tr(e- Li (- p(1)h ^ p(2)h2 ) • f[g(\) D„ X(A) D^ 1 ) , (2.16) 

^ n=l ' 



where 3x3 diagonal matrix D n is given by 

D n = exp h 3 + ^ 2) h 2 Y 



A matrix is (E^)^ = 6 i>m S j>n . 



h a = E aa --t. 



For our purpose to construct the Baxter equation for the discrete Boussinesq equation, we use 
another transfer matrix which is associated with a representation 3* . Generally the monodromy matrix 
for this type of representation can be constructed from one for the vector representation by the fusion 
method [24], and in the s£^ case the Lax matrix is simply defined from L n (A) oc nL n (A)) ^ , 
where t denotes a transpose [|25|] ; 



L n (A) = g(X) 



\ 



^A0« + lA^ 



/e A -1 0\ 



e A - 
\-l e A J 



(2.17) 



5(A) • C - 1 • X(A), 



where as before 77(A) is a normalization function. One sees that each element of C n _1 corresponds to 
weights of 3*. Correspondingly the monodromy matrix and the transfer matrix for this representation 
is respectively defined by 



T(A) = ;QMA), 



t 2 (A) = TrT(A). 



(2.18) 



We note that under the boundary condition (2.15) the transfer matrix i 2 (A) has a form, 



t 2 (A) =Tr(e Li (- pWh ^ P(2)h >) • f[g(\) e** D^ 1 X(A) D n ) 



(2.19) 



The commutativity of this transfer matrix follows from eq. (2.12) as follows. We have the funda 



mental commutation relation for T(A) from eq. ( |2.12| ) as 



._i\tit 2 J_ JL A J_ /. ,_i\tit 2 
R(A -/i)-) ) T(A)T(/i) = T( M )T(A)((R(A-/x)) 



(2.20) 
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where t a denotes a transpose over the a-th space. We see the commutativity of the transfer matrix; 

MX) , t 2 (n)\ = 0. (2.21) 



Furthermore eq. (2.12) also gives 



l 



(R(A - /,)) t2 ) • T(A) T( M ) = T( M ) T(A) • ((R(A - /,)) t2 
which proves the commutativity among two transfer matrices, 

MX) , t 2 (/i)] = 0. 



(2.22) 



(2.23) 



2.2 Fundamental £ Operator 

We shall construct the generating function of the local integrals of motion for the affine Toda lattice 
field theory. 

For the Heisenberg operators p and q satisfying 

[p, q]= -2i7, 
we define the fundamental C operator as [^2|]p| 

£(A;g,p) 



e 7 (g) e 7 (p) 



(2.24) 



$ 7 (A + p) $ 7 (2A + g) $ T (A-g) $ 7 (A-p)' 
Here A is a spectral parameter, and commute with every quantum operators, and the functions & y ((p) 
and © 7 (</?) are respectively defined by 



&~f(<p) = exp 



-upx 



dx 



lu+io 4 sh(jx) sh(vrx) x 
@ 7 (ip) = $ 7 0) ■ $ 7 (-p). 



(2.25) 



(2.26) 



The integral «$ 7 (v?) was first introduced by Faddeev [27] as a non-compact version of the quantum 
dilogarifhm function. See Appendix |A] for properties of these functions, and we especially recall that 



the function 7 (</?) is the Gaussian ( |A.3[ ). By use of the difference equation ( |A.4[ ), we have y22| ] 

1 



£(A; g - 2 i 7 , p) = £( A; g, p) • (e A + e"*-^ + e"^ 



1 



e -A + e A-i 7 +g + e i 7 +p 



) e A-i 7 +g _)_ e -A-i 7 +p 

(1 + e -A-i7+9 + e A+i 7 +P) . £( A; ^ 



(2.27) 



£(A; g, p - 2 i 7 ) = £(A; g, p) • (1 + e A+i ^ + e" 



A-i 7 +p) 



1 



1 -)_ e -A-i 7 +g _|_ e A-i 7 +p 



•(e A + 



-i 7 +<? 



+ C 



-A— i 7 +p 



) • £(A;g,p). 



(2.28) 



" Strictly speaking, the operator ( |2.24| ) is a non compact ana logue of that was introduced in Ref. [22J, i.e., we replace the 
g-exponential function (A. 1 ) with the Faddeev integral integral ( |2.25[ ) J^], 
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As two operators Acfffl and A<pn^ are the Heisenberg operators, [A</>.„^ , A(pffl] = 2ij, we 
denote for brevity 

C n (X; <j>) = £(A; A0« , A^ 1} + A#). (2.29) 

We collect useful properties of this operator in Appendix |b[ This operator acts as the intertwining 
operator for the vertex operators satisfying [ |2l| ] 

(e^ + e A e*»+i) • C n (X; <j>) = C n (X; <j>) ■ (e A + e*&) , (2.30) 

for a = 0, 1, 2. The duality ( A2) of the integral ( |2.25| ) indicates that the operator C n (X; <ft) also 
intertwines the dual vertex operator; 



lS a ) 7T \ TV jX a ) 

e7*» +e^e^) ■ C n {X;<j)) = £ n (A;< 



e^" + e^™ +1 )- (2.31) 



What is important here is that the fundamental C operator Q2.29D itself satisfies the Yang-Baxter 
equation [[28]], 

£ n (A; 4>) £ n +i(A + /u; 0) £„(//; 0) = £ n+ i(//; £ n (A + /x; 0) £„+i(A; 0). (2.32) 
Thus when we set the Q operator as 



Q(\)=Hf(\)-C n (\;, 



(2.33) 



with a normalization /(A), we see that they commute between themselves with arbitrary spectral 
parameters; 

[2(A) , Q(jj)] = 0. (2.34) 
Furthermore the fundamental C operator also satisfies j j22| ] 

£„(A - fi; 4>) L n+ i(A) L n (/j,) = L n+ i(/j,) L n (A) C n (X - 4>), (2.35) 



which proves that the transfer matrices t\(X) and £2 (A) commute with the Q operator ( |2.33[ ), 

[2(A) , hQi)] = [2(A) , * 2 (m)] = 0. (2.36) 

In conclusion we have a commutative family; the transfer matrices ti(X) and £2 (A), and the operator 
Q(A). We note that intertwining relations ( 2.3Q| ) and (2.31) proves the commutativity between the Q 
operator and the screening charges, 

[2(A) , Q {a) ] = [2(A) , QW] = 0. (2.37) 

As a consequence, the Q operator is a generating function of the local integrals of motion of the affme 
sis Toda field theory [[21]], whose Hamiltonian and its dual are 

2 2 



a=0 



a=0 



See Refs. [29|,|30| for the geometrical interpretation of the screening charges in the classical limit of the 
sine-Gordon system. 



8 



3 Baxter Equation 



We shall show that the Q operator defined in eq. (2.33) plays a role of the Baxter Q operator for the 



transfer matrices (2.11) and (2.18) of the quantum discrete Boussinesq equation. 



Step We prepare some formulae of the fundamental C operator (2.29); using elements of matrix 



C n = diag(cn ■* , c£' , ) in eq. ( ^9| ), we have 

2 



e- A -^-£ n (A + -i 7 ;< 



4 2) (A;<A)-(c« +C ( 2 )e A + c( 3 )e- A 



e-3^(e 3A -l)-£ n (A-^i7;, 



,(2) 



4 1 )(A;^)-e A -4 2 )(A 



cl 3) ) -(4 2) (A;0)-e A -4 3 )(A 
cW) _1 -(4 3) (A;0)-e A -£«(A 



Here we have defined £„ (A; 0) by 



which, using the commutation relations (|2.2|), are explicitly written as 

(A; 4>) = £(A; A0« - | i 7 , A0« + A<^ 2 ) + H i 7 ), 

4 2 ) (A; 0) = £(A; A0« - | i 7 , A0« + A# - | i 7 ) , 
£f (A; <f>) = £(A; A0« + ^ i 7 , A0« + A# + | i 7 ) . 



See Appendix ^] for proof of eqs. (3.1) and (3.2) 



(3.1) 



(3.2) 



Step 1 We first consider the product of the Q operator ( |2.33| ) and the transfer matrix ti(A). As we 



know from the commutation relations (^2f) of the lattice free field that a difference field A(j)^ does 
not commute only with nearest neighbor A^j^, we get 

2(A) • h(X) = Tr J] /(A) 5 (A) £ n (A; 0) C n+1 X(A) 

n 

= Tr J] /(A) g(\) M" 1 C„-\ £ n (A; 0) C n+1 X(A) C n M. (3.3) 
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Here a matrix M denotes a gauge-transformation. We substitute 

M = 1 1 

U h 



into eq. (p.3|), and compute each element directly. We see that both (2,1) and (3,1) components van- 



ishes due to eqs. < J3- 1[ ) - Thus we find that a trace of the product of 3 x 3 matrix reduces to a sum of the 
product of (1,1) component and a trace of the product of remaining 2x2 matrices; 

Q(A) • h(X) = f[ /(A) g(X) ■ C n (X + | i 7 ; <f>) + Tr f[ /(A) g(\) A n (A), (3.4) 

n n 

where 



A n (A) 



-n i; (A;0)-e" A + £^(A;</))-e A — (A; 0) + £„ (A; 0) • e \ /c}i 
-^(A^-e^ + ^W) -£^(A; 0) + £l 3) (A; 0) • e A j I cl 3) 



We now rewrite above 2x2 matrix A n (A). Using eqs. (3.2), we have 
A„(A) = e-^(e 3A -l) I 



r W e -A .(1) J " v 3 " r/ \ _(3) 



e" A ( e 3A - 1) 



'd 3) (A-|i 7 ^ 



A lj (A-|i 7;i 



In the last equality, we have used the commutation relations (2.2) and a condition ch^ c4 c^ 3 ^ = es 17 . 
As a result we get 

Q(A) • ti(A) - fj /(A) 5(A) e~ A -^ • C n {\ + \\ T A) 

n 

X^y' + ^y 1 -^ -(civ 



Ti-n/(A)5(A)e- A (e 3A -l) 



£l 3) (A-|i 7 ^) 

^(A-fiT^) 



(3.5) 



We keep aside this identity for a while. 
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Step 2 In the same manner with Step 1, we study a product of the Q operator and the transfer matrix 
i 2 (A). We have 

Q(A)-* 2 (A) 

= Tr /(A) 5(A) £n(A; 0) C - + \ X(A) 

n 

= Tr H /(A) 5(A) C n+1 £ n (A; 0) C - + \ X(A) C - 1 

n 

= TrJ]/(A)5(A)M- 

V-^V'-^W) 
In the last expression, we have introduced matrix M for gauge-transformation. When we substitute 

'I 



\ 

ri 3) )-V-£l 2) (A;^ 



M _1 . 



M 



1 1 1, 



we see that both (3,1) and (3,2) elements vanish due to eqs. (3.2), and that a trace of the product of 
3x3 matrices can be rewritten as a sum of the product of (3,3)-element and a trace of the product of 
2x2 matrices; 

2(A) • h(X) = J] /(A) g(X) e~& (e 3A - 1) ■ £ n (A - - i 7 ; 0) + Tr J] /(A) 5(A) B n (A), (3.6) 



B„(A) 



-(^^•^(A;, 



where 2x2 matrix B„(A) is given by 

(c^V 1 •4? ) (A;£) rfV^-^t^) + ri 3 ))- 1 -4 3) (a ;i 

We see that by applying eqs. ( |3~^ ) to (2,2)-element of above matrix B n (A) we have 



TrJJ/(A)5(A)B n (A) 

n 

'(ci 1) )" 1 e A • d 3) (A; <f>) -(c { n ] r l • d 2) (A; 0) 

(cf)- 1 -^^) ((^VV + ^V 1 ) -4 2) (A;0 
where the invertible matrix M is included again for gauge-transformation. Substituting 



T*n/(A)5(A)M. 



•M 



M 



1 r 
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and erasing C n (X;<p) using eqs. < p?2\ ), we get 
Tr[]/(A)5(A)B n (A) 

n 

\ (cf)- 1 (c^V'eV V Cn (X; <j>) t 

As a result we obtain an identity, 



2(A) • t 2 (A) - f[ /(A) 5(A) e"^ (e 3A - 1) ■ ^(A - | i 



7; 



Trn/(A) 5 (A) ( , • J - M 



(Cn 



Step 3 By comparing eq. (3.5) and eq. (3.7), we find that 2x2 matrices in the right hand sides have 
same form with each other. By erasing these 2x2 matrices, we obtain an identity for operators t\{X), 
f 2 (A), and Q(A); 

QW ■ tl(A) " (5 ?<rrfw) s(x) e "~ il7 ) • Q(x + 1 i7) 

(fj /(A)g(A) _ A/3A A 2 - 1( ,, 2 - i 

= (n /(A _ |i7)5(A _ |i7) e < e -l)J«A-s>7)-«A- i .7) 

- (n7^|^9A)e- A ' l ' 7 (^-l)(e 3A - 2i ''-l)) <2(A-|i7)- 

After setting normalization functions as 

/(A) = 1, 
5 (A) = e A +^, 
5(A) = e^ (e 3A+2i ^-l), 
we see that the operator Q(X) solves the difference equation, 

Q(A + 2i 7 ) - ti(A + - 17) • S(A + 1 17) 

+ * 2 (A + 1 17) • Q(A + I i 7 ) - A(A + | i 7 ) • Q(A) = 0, (3.8) 
where the function A (A) is defined by 

A(A) = ((e 3A -l)(e 3A+2i ^-l)) L . 
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As we know from eq. (2.36) that the transfer matrices ii,2(A) and the Q operator commute each other, 



the operator valued difference equation (3.8) can be regarded as the third order difference equation in 
usual sense, once we apply both hand sides to simultaneous eigenfunction. 



We expect in general the third order difference equation (3.8) has three linearly independent solu- 
tions, which we write Q a (X) for a = +, 0, — . We then introduce the function V(X) as the quantum 
Wronskian of solutions; 

§*■ 

2 • , 



v ab {\) 



Qa(\) Qa(A-§i 7 ) 



(3.9) 



QbW Q 6 (A-§i 7 ) 

for a, b = 0, ± (we thus have three independent functions V a b (A))- It is straightforward to see that the 
function V a b(X) satisfies the third order difference equation, 



P(X + 2i 7 ) - t 2 (X + | i 7 ) • V(X + | i 7 ) 

+ A(A + ^i 7 )t 1 (A + ^i 7 )-P(A + ^i 7 )-A(A)A(A + ^i 7 )-P(A) 



0. (3.10) 



From two difference equations (3.8) and (3.10), we get the Baxter equation for the discrete Boussi- 
nesq equation; 



h(X) = A(A 



2. V(X 
3 17J V{X) 



|i 7 ) , P(A + |i 7 ) G(A-|i 7 ) , S(A+|i 7 ) 



+ 



h(x) 



P(A + fi 7 ) 



+ A(A) 



V(X) Q(X) 
V{X) S(A + |i 7 ) 



+ 



2(A) 



P(A + |i 7 ) 



2(A) 



+ A(A) 



2(A-|i 7 ) 
2(A) 



Once we have a set of functional equations ( p. 1 1| ) we can follow a strategy in Ref. 
spectrum of these transfer matrices; we suppose 



Q(A) = nsh(A-A( 1 



V(X) = Y[ s h(X-X^), 



(3.11a) 

(3.11b) 
to obtain the 

(3.12) 



and substitute them into eqs. (3.11). From the condition of analyticity of the transfer matrices t a (X) 
we get the nested Bethe ansatz equations. The precise analysis of these equations is for future studies. 
As a consequence, we have diagonalized the transfer matrices ti,2(A) by use of the auxiliary functions 
Q(A) and V(X). 



4 Duality 



We have seen that the Q operator (2.33) satisfies the third order difference equation (3.8) whose 
coefficients are given by two transfer matrices ti(X) and £2 (A) for the quantum discrete Boussinesq 
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equation. As was noticed in Appendix [A| the Q operator has a property of duality which follows from 



eq. (2.31) and eq. (A2); 



7T- 



1 



7 



A — >• — A. 



(4.1) 



7 7 7 

Above duality implies the existence of the dual Baxter equation; we have the dual third order differ- 
ence equation, 

Q(A + 2ivr) - ui(A + - Itt) • Q(A + ^ ivr) 

+ u 2 (A + | ivr) • Q(A + | ivr) - A(A + | ivr) • 2(A) = 0, (4.2) 



where 



A(A) 



e 3 T A — 1) (e~ 



(3A+2i7r) 



The operator Q(A) is as same with eq. (2.33), and the transfer matrix u a (X) is defined by 
ui(A) = TrU(A), U(A) = J]M n (A), 



n 2 (A) =TrU(A), 
with the dual Lax matrix M(A) and M(A), 



U(A) = []M n (A), 



M n (A) = h(X) 



C" 



V 



/ P ^(-|A^ + iA^)) 



C 



^(-|Atf)-lA^ 2) ) 



M n (A) = h(X) 



g7 V 3 



|A^)-|A^) 



i? (|A^) + iA^; 



(4.3) 



(4.4) 



/ -A — -A 
i -A 

1 eT 
e t 1 



/ ~ A 1 



eT 



-1 



To get the dual Baxter equation, another operator V should be slightly different from eq. (p.9|), 



V ab (\) 



Qa(A) Qa(A — | 171") 

Q 6 (A) Q 6 (A-|i7r) 



(4.5) 



where Q a (A) denotes three linearly independent solutions of the difference equation (12). This solves 
the third order difference equation, 

P(A + 2ivr) - u 2 (A + | i vr) ■V{\ + \\n) 

+ A(A + \ ivr) m(A + \ ivr) • P(A + \ ivr) - A(A) A(A + \ in) ■ V(X) = 0. (4.6) 

o o o o 
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As a result, we obtain the dual Baxter equations, 



V ' v 3 ; v(\) V{\) 2(A) Q(A) 



. 2 (A) = £^±4^ + A(A)^ 



n\) 2(A + fi7r) Q(A-ji7r) 



2(A) 



+ A(A) 



2(A) 



(4.7a) 



(4.7b) 



P(A + f ivr) v ' P(A + |ivr) 

from which we get the nested Bethe ansatz equations only from the analytic property of the transfer 
matrices u a (X). 

We remark that this duality of the Baxter equation also appeared in studies of the discrete KdV 



equation [31J, where given was the interpretation from the viewpoint of the algebraic geometry. 



5 Continuum Limit 



We briefly study the continuum limit of the quantum discrete Boussinesq equation. In a continuum 
limit, the lattice free field (p^ is replaced by (fi^ (x), and the commutation relations ( jL2| ) then reduce 
into 



\^ a \x), ^ b \y)] = C ab i 7S gn(x-y), 
where C is the Cartan matrix of sis, 

2 -1 
-1 2 



The quasi-periodic boundary condition ( |2.15| ) denotes 

<f>W(x + L) = (f) {a \x) + L-iP( a \ 



which indicates the mode expansion of the free fields as [32] 



(b) 



4>^(x) = iQ^+iP^x + y ^e 1 



(5.1) 



We can take the continuum limit of the transfer matrices ii,2(A) from eq. ( |2.16| ) and eq. ( |2.19[ ). 
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After a gauge-transformation, we get 

{g(X)e x+1 ^y L -t 1 (X) 



Tr 



^(pWhi-pWhs) p 



exp 



dx 



e - 2A- T i7 e-* (1) ( x ) e - A+ T i7 e* (0) ( x ) \ 

-A+|i7 e (1) W e -2A-ii7 e -^ 2 >W 





(5.2) 



(5(A)e A+ ^)- L -t 2 (A) 



Tr 



3 iL(-p(i)hi+pC2)h 3 ) p 



cxp 



-A-li 7 



d.r 



/ 













(5.3) 



where & denotes the path operator ordering. One finds that the transfer matrix £2 (A) in eq. (5.3) 



essentially coincides with that was proposed in 32 as a quantized Drinfeld-Sokolov reduction. 



Correspondingly a difference of the free fields (2.4) becomes the differential of the free field, 





A6^ 



dx 



and the Baxter Q operator for the quantum Boussinesq equation can be given by the path ordering of 



the operator C n (X; </)). As a result we recover the third order difference equation equation (3.8). At 
this stage we are not sure how to relate the Q operator with the universal 1Z matrix. 



6 Generalization: TV-reduced Discrete KP Equation 



As the Boussinesq equation ( fT"2| ) is the 3-reduced KP equation, we can introduce the quantum N- 
reduced discrete KP equation which is related with the affine s£n Toda field theory. In the classical 



case, the Lax matrix for the discrete A r -reduced KP equation is given by [33, 34] 



C(x) 



1 





V 



JV-2 



VV n+N-2 X 



1 






N-3 



w. 



(3) 








n+N-3 



{-l) N - 2 xW { n N + - l) (-l)K-i W { n N) \ 








(6.1) 
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where the generators of the lattice Wjy algebra are defined as 



= — ^ ( E ex P (xfii-i + xSi- a + • • • + xfe 13 



+ £ exp( X ^_ 1 + X ^_2 + --- + X^ ) ) , for S = 2,3,...,iV-l, 

l<il<i a <— <ia<]V-l / 



exp 



'JV-1 N 
v i=l > 



W-l / JV-1 N 

II l+E eMx^l) 

k=0 \ i=l / 



(6.2) 



[4 a) , 4 a+1)] 



Here Xrf 3 is defined in terms of the free fields <j>$ (for a = 1, 2, . . . , N — 1) as 

6=1 

and the quantum algebra for the lattice free fields reads as 

[^),#] = 2i 7 , for m>n, 

J i7, for n < m, (6-3) 
1 —17, for n > m. 

We note that the transformation ( |6.2[ ) is the so-called Miura transformation. We set 

JV-l 

# = -E^ 

a=l 

and the screening charges are defined by (for a = 0,1, . . . , N — 1) 

Q (a) = ^ e ^ ; QW=2e^ ) . (6.4) 

n n 

The sum of the screening charges, J2a=o Q^> corresponds to the Hamiltonian of the affine Toda 
field theory. 



In terms of the free fields, the quantum Lax matrix is then defined by [ p4| 



L n (A) = exp^E Kx ( n ] ) • X(A), (6.5) 



where N x N matrices h a and X(A) are defined by 

JV-l N-l 

K = E , a - - 1, X(A) = E ^ C fc , c = Ei,^ + E E *+^" 

fc=0 fc=l 
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The spectral parameter A is related with x in eq. ( |0| ) by 

As usual the monodromy and the transfer matrix are respectively given by 

T(A) = J]L n (A), h(X) = TrT(A). (6.6) 

This transfer matrix generates the integrals of motion of the discrete N -reduced KP equation. 
The fundamental £ operator for sijq case is defined by [ ]22| ] 

1 1 1 



£n(A; ' 



<MA + X ^ -1) ) $ 7 (2A + X ^~ 2) ) $ 7 ((iV-l)A + x^ 



x e 7 ( x ^) Qy(x ( n ] ) e 7 (xf" 1} ) 

^(A-x^) $ 7 (A-^ 2) ) $ 7 (A- X ^ 1} )' 
This operator not only satisfies the Yang-Baxter equation ( [2.32D , but fulfills the intertwining relations 



for the (dual) vertex operators Q35Q ; 



(e^ + e A e*&) • £ n (A; 0) = £ n (A; 0) • (e A + e*&) , (6.8) 
fe^ n +e7 A e7^+ij .£ n (A ;( /.) = £ n (A;(/))- fe7 A e 7^ +e ?^+iJ, (6.9) 
for a = 0, 1, . . . , N — 1. We can easily see that the Q operator 

Q(\) = l[£ n (\;0), (6.10) 

n 

commute with the screening charges 

[Q(A), Q^] = [Q(\), QW]=0, (6.11) 

which proves that the Q operator generates the local integrals of motion for the discrete analogue of 
the affine s£n Toda field theory. 

As we have clarified for sis case in preceding sections, the generating function of the local inte- 
grals of motion for the affme Toda field theory becomes the Q operator for the quantum Boussinesq 



equation. Thus it is natural to conjecture that the Q operator (6.10) satisfies the A^-th order difference 
equation, 

N-l 



N — k N — k 2 

Q(A + 2i 7 ) - ^ A + 2 —5^7) • 2( A + 2 — aT" i7) + A(A + N i7) " Q(A) = °' 



N " v A r N 

k=l 



(6.12) 
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with some function A (A). Here the transfer matrix ii(A) is as eq. ( |6.6D and corresponds to the vector 
representation. Other transfer matrices ifc>i(A) are for [l fc ], and can be constructed by the fusion 
method from the transfer matrix ti(A) with a suitable normalization. Especially we have 

/ 2V-1 \ 

t N -xW oc Tr II ex p(- E • X ( A ) ' 

\ n a=l / 

with X(A) = 1 — e~ A C" 1 . The duality of the Q operator also suggests the dual difference equation, 

1 lz} N — k N — k 2 

Q( A + 2ivr) - + 2 — ^— ivr) • Q(A + 2 — — ivr) + A(A + — Itt) • 2(A) = 0, 

k=l 

(6.13) 

where the transfer matrices itfc(A) follows from ifc(A) under a duality transform (fO|). 



7 Concluding Remarks 



We have explicitly constructed the Baxter Q operator for the quantum discrete Boussinesq equation 
which has U q {st^). As far as the author knows, the explicit form of the Q operator has been known 
only for the quantum integrable systems associated to si2, and our Q operator is the first example 
which solves the third order difference equation. In our construction, the Q operator originates from 
the generating function of the local integrals of motion for a discrete analogue of the affine s£s Toda 
field theory, and we hope this can be applied for other Lie algebras, to which some attempts have been 
made [f36|]. Important structure is that it has a kind of duality, and that the Q operator satisfies the 



dual Baxter equation. It was shown [31] in the case of the discrete KdV equation that the dual Baxter 
equations can be interpreted as the duality between homologies and cohomologies of quantized affine 
hyper-elliptic Jacobian from the view point of the algebraic geometry. It will be interesting to give 
such interpretation in our case. In such studies the separation of variables (see Ref. [37| for review, and 



Ref. 38 for st^ case) will be useful. 



Also in sl2 case, the universal procedure to derive the Baxter equation was proposed [14], where 



the Q operator was constructed from the universal 1Z matrix with the infinite-dimensional (/-oscillator 



representation. This may be true for s£n case, but we are not sure how to relate our Q operator (6. 1C) 
with the universal 1Z matrix. 

The Baxter equation for sl2 appeared in studies of the second order ordinary differential equa- 
tions; the spectral determinants for the Schrodinger equation is identified with the Q operator for the 
CFT [^9|-4TTj]. Such correspondence is also studied for the higher order differential equation [[l||, 43]. 



This "ODE/IM correspondence" seems to come from a similarity between the Stokes multiplier and 
the Lax matrix, and our results will help further investigations on these topics. 



19 



Acknowledgement 

The author would like to thank Rei Inoue for stimulating discussions. 



20 



A Quantum Dilogarithm Function 



We review a q-deformation of the dilogarithm function ( 2.25| ); 



&-y((p) = exp 



-llfX 



d.r 



„ i0 4 sh(jx) sh(7rx) x 



(2.25) 



This integral was introduced by Faddeev []27|], and it conesponds to a non-compact analogue of the 
(/-exponential function, 



oo / oo 

iV x T9 2n+1 iu) = exp 

n=0 



S q (w) = H{l + 



(-l) k w k 



(A.l) 



We list properties of the integral ( 2.25 ) in order [ 44 ] . 



Duality 



(<p) = * 7 (-P)> 

V 7T 



(A.2) 



• Zero points, 

zeros of ($ 7 (^)) ±1 = j=Fi((2 m + 1) 7 + (2 n + 1) tt) 

• Inversion relation, 

& 7 (<p) = $ 7 (y>) • * 7 (-¥?) 



m,n £ Z> 



(2.26) 



exp 



1 ^ 2 7T 2 + 7^ 



(A.3) 

2i 7 v 2 6 yy 

A reason why we use notation "0" is that it is proportional to the Jacobi theta function when 
we replace the integral 3> 7 (<£>) with S q (w), S q (w) S q (w~ 1 ). 



• Difference equations, 



$ 7 (p + i 7 ) 1 



$ 7 ((/?-i7) ~ 1 + e^' 
The duality ( |A.2[ ) gives another type of difference equation, 

$ 7 (y? + i7r) _ _J_ 
$ 7 (v?-i7r) ~ i + e ^' 



Pentagon relation [ }44] , |45p 

where j5 and q are the Heisenberg operators satisfying a commutation relation, 

[p , q] = -217. 



(A.4) 



(A.5) 



(A.6) 



(A.7) 



This identity is a key to apply the integral ( [2.25 ) to the quantization of the Teichmiiller 



space Q45|] and to the construction of the invariant of 3-manifold [46, 47]. 
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B Properties of the Fundamental C Operator 



We study the property of the fundamental C operator (2.29). 



B.l Z 3 Symmetry 



By construction the free fields 4>n \ 4>n \ and 4$ 



i>n — fin correspond to roots a±, e*2, and 



ocq. As the fundamental C operator C n {\; <p) is defined as the intertwining operators for these Z3 



symmetric vertex operators (2.30), we can suppose the C operator itself has the Z3 symmetry. This 



can be checked by recursive uses of the pentagon identity ( |A.6| ) as 

£ n (A; d>) = ($ 7 (A - A<j>&>) $ 7 (A - A0«) $ 7 (2 A + A^ 1 )) $ 7 (A + A0« + A#)' 
xG 7 (a4 2 ))G 7 (A4 1 )) 
\(X + A^ 1 ) + A^ 2 )) <D 7 (A - A#) d> 7 (2 A + A#) <D 7 (A - A^ 1 )))^ 

xG 7 (A4 1 )+A4 2 ))e 7 (A^ 2 )) 
\(X - A^ 1 )) $ 7 (A + A0« + A#) * 7 (2 A - A^ 1 ) - A#) * 7 (A - A^ 
xG 7 (A4 1 ))e 7 (A4 1 ) + A4 2 )). 

As a result, we see that 

£ n (A; 0) = £(A; A0«, -A<^ 0) ) = £(A; A</>^,-A<f>^) = £(A; A<^ 0) , -A#). (B.l) 



B.2 Unitarity 



Recalling the fact that the Q function (A. 3) is proportional to the Gaussian, we can check that the 
fundamental C operator satisfies the unitarity condition, 



£ n (\;4>) • £ n (-A;< 



e 2i ~< 



^- A 2 



(B.2) 



C Proof of Eqs. (£Tp - 



We sketch the outline of a proof of eqs. (3. 1) - (3.2) 
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To prove eqs. (JO|), we first compute as follows; 
1 



£ n (A + §17;. 



. £ (i) (A . 0) = ?7(A + |i7 - - A^ 2) ) 7 (A<^ - |i 7 ) 



9 7 (A^ 1} +A</>. 



e 7 (A^ 1} ; 



3> 7 (2 A + A^ + | i 7) 9 7 ( A^ + A<^ 2) + § i 7) 



$ 7 (2 A + A<$> - § i 7 ) $ 7 (A - § i 7 - A^ 1} - A<^ ; ) 



$ 7 (A + |i 7 - A^ 



(2K 

n J 



e 7 (A^ 1} + A^) 

e 7 (A0« + + |i 7 ) _ i A ^)-i i7 



(2), 



2A+ii 7 +A0L 1) 



1 + e 



$ 7 (A-|i7-A^ 1} -A^ 2) ) 



A+|i7 + e 2A+|i7+Atf ) + R A^ 1) +A^ 2) +i 7 



8 7 (Atf > + A# + |i 7 ) e i A ^)-|i 7 



e 7 (A^ 1} +A^ + 2i7) 



(2) 



C3 



e A+ii 7 + g2A+|i 7 +A0„ ' + e A^+A^ 



(CI) 



This proves the first equality in eqs. (|3JJ). Remaining equalities follow from eqs. ( |2.2V| ). 



£n(A 



To prove eqs. (|3.2|), we shift parameters in eq. (|3JJ); 

1 a " • ,m |i 7 ,A^) + A^)-| 



£(A; A0« + § i 7 , + A# - \ i 7 ) = e~fr (c« e 2 ^ + + e> 



(C.2) 



Using eqs. (3.1) and above identity, we can prove the first equality in eqs. (3.2) as follows; 

£«(A; <t>) ■ - C%\\; 0) = 4 2 )(A; 0) • ((c« + c£V + c^e^) • 



r W e \ , J 2 ) p -A , r (3) 



4 2) (A;0)-(< 



2A 



e- A ) • a 



(e 3 ^ - 1) • c( 2 ) • £(A; A# + | i 7 , A0« + A# - | i 7 ) 
1 



J!)p2A 1 J 2 ) 1 J 3 ) p A 



( e 3A_ 1)e -|i 7 . c (2). £n(A __ i7 .^. 



Other equalities directly follow from the Z3 symmetry (B.l). 
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